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THE ORIENTATION OF THE NUCLEAR OBSCURER OF THE AGNS 
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ABSTRACT 

We examine the distribution of axis ratios of a large sample of disk galaxies hosting type 2 AGNs 
selected from the Sloan Digital Sky Survey and compare it with a well-defined control sample of non- 
active galaxies. We find them significantly different, where the type 2 AGNs show both an excess of 
edge-on objects and deficit of round objects. This systematical bias can not be explained by a nuclear 
obscurer oriented randomly with respect to the stellar disk. However, a nuclear obscurer coplanar 
with the stellar disk also does not fit the data very well. By assuming that the nuclear obscurer having 
an opening angle of ~ 60 degree, we find the observed axis ratio distribution can be nicely reproduced 
by a mean tilt angle of ~ 30 degree between the nuclear obscurer and the stellar disk. 
Subject headings: galaxies: statistics — galaxies: Seyfert — galaxies: nuclei — galaxies: spiral 



1. INTRODUCTION 

In the unified AGN model, a type 2 AGN is seen from 
an angle close to the plane of the central accretion disk. 
On this plane, the accretion disk and the broad line re- 
gions are obscured by an outer molecular torus, while a 
bi-polar jet em anating from the nucleus is perpendicu- 
lar to the plane (|AntonuccilH993l ). The type 1 AGNs are 
then these objects that we are able to see directly into 
the central regions. 

Several studies have revealed that the angle between 
the orientation of the kiloparsec-scale jet and the nor- 
mal to the host galaxy plane has a wide/random 
distribution based on the samp l es of dozens of ra- 
dio galaxiesdSchmitt et all 119971; iKinnev et all 120001: 
IGallimore et al.l 120061) . These results lead to a conclu- 
sion that the internal gas fuelling the black hole may not 
be co-aligned with the outer gas forming the stellar disk. 

On the other hand, the AGNs are known to be biased 
sample of inclinations of stellar disks for a long time. 
Both the optically selected (mostly type 1) and the soft 
X-ray selected Seyfert galaxies show a bias against the 
edge-on galaxies, whi le the hard X-ray selected sample do 
not show such a bi as(IKeellll980t lLawrence fe Elvislfl 982; 
ISimcoe et aHfl997t [Zhang et al.ll2009l) . IMcLeod fe Riekd 
(119951) found an excess of face-on galaxies hosting 
both type 1 and type 2 AGNs and explained it with 
a substantial amount of ob scuring material coplana r 
with the stellar d i skfsee also Ide Zotti fc Gaskeli 119851 ) . 
IMaiolino fc Riekel (|1995l ) showed that the intermediate 
type AGNs(1.8,1.9) are mostly found in edge-on galax- 
ies and concluded a 100 pc-scale torus coplanar with the 
stellar disk. 

In this study, we use a large type 2 AGN sample se- 
lected from the Sloan Digital Sky Survey (jYork et al.l 
120001 SDSS) to check whether the stellar disks of their 
hosts are biased to edge-on views. We then further check 
whether this bias is better explained by a dust layer 
coplanar with the galaxy plane or by an intrinsic align- 
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ment between the nuclear obscurer and stellar disk. This 
paper is organized as follows. In section 2, we study the 
axis ratio distribution of the AGNs and make compar- 
isons to that of a well-defined sample of non-active disk 
galaxies. With the distributions of the shape paraments 
of the spiral disks concluded in Section 3, we then build 
toy models to quantify the biased axis ratios of the AGN 
hosts in Section 4. Finally, we make brief conclusions 
and discussions in Section 5. 

2. THE INCLINATIONS OF THE STELLAR DISKS 
OF THE TYPE 2 AGN HOSTS 

Our type 2 AGNs are selected from the com- 
plete spectrosco pic main galaxy samp le of the data 
release 7(DR7, lAbazaiian eTail 12009 ) of SDSS us- 
ing the empirical relat ion between th e emis sion 
line ratios proposed by iKauffmann et al.l (|2003l )(see 
their equation 1). The emission line data are 

take from th e MPA-JHU DR7 release of spectrum 

measurements( http : / / www . mpa-garching . mpg . de/ SDSS/DR7/) We 
use the SDSS parameter ' fracDeV 1 and criteria 
fracDev < 0.5 0, which ensures that the galaxy flux is 
dominated by an exponential(disk) component, to select 
the AGNs hosted by disk galaxies. The number of the 
type 2 AGNs hosted by disk galaxies is then 32,618. 

The inclination of a stellar disk is informed by its ap- 
parent axis ratio b/a. If a disk is circular and thin, then 
b/a = cosi, where i is the inclination of the disk (the an- 
gle between the line-of-sight and the normal of the disk). 
If this disk is viewed from arbitrary directions, the re- 
sulted b/a follows a uniform distribution between and 
1. 

For a more realistic triaxial ellipsoid model, the ob- 
served b/a is further dependent on two shape param- 
eters, the disk height 7 = C/A and the disk elliptic- 
ity e = (1 — B/A), where A,B,C are the major, mid- 



1 fracDeV is defined in Abazaiian et al. (2004) through Fcom- 
posite = fracDeV FdeV + (1 - fracDeV) Fexp, where Fcompos- 
ite, FdeV and Fexp are the composite, de Vaucoulcurs and expo- 
nential fluxes of the object respectively. 
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die a nd minor axis of the ellipsoid respectivelv (|Binnevl 
1985). The b/a distribution is then dependent on both 
of the shape parameter (7, e) and viewing angle distri- 
butions. The galaxies with different physical properties 
(e.g. mass and size) are shown to have differ ent intrin- 
sic shape parameters ()Padilla fc Strauss! 12008'). In addi- 
tion, the sample selection effects could also introduce bi- 
ases into the measured b/a distribution of a galaxy sam- 
ple. For example, the seeing would make the galaxie s 
with small apparent sizes look rounder (jShao et alJl2007[ ). 
Therefore, to quantify the inclinations of the disks of the 
AGN hosts, their shape parameters must be well-defined 
or controlled. 

The AGN host galaxies are known to have distinc- 
tive physical properties. F or example, the AGNs r eside 
mostly in massive galaxies (Kauffma nn et al.l I2003T ) and 
the strength of th e nucleus activity is also correl ated with 
its host propertv (|Kauffmann fc Heckmanll2009H . For our 
sample of AGNs hosted by disk galaxies, we find that 
they are biased towards these hosts with large bulge 
component, high concentration and old stellar popula- 
tion. To quantify the sample properties of these AGN 
hosts, we build a control sample of galaxies (without 
AGN phenomena) from the same galaxy catalogue of 
SDSS, which are selected to have the same sample size 
and the same distributions of stellar mass, redshift, size, 
concentration, fracDev and D n (4000) (the 4000A break 
spectra index, a rough stellar age indicator) as the AGN 
hosts. We believe that this control sample of non-active 
galaxies shares the same stellar properties and selection 
effects(e.g. redshift distribution) as the AGN hosts so 
that they would have the same b/a distribution if their 
viewing angles were both random. 

We compare the b/a distributions of the AGN hosts 
to the control non-active galaxies. The result is shown 
in the top panel of Fig. [TJ The AGN sample (solid his- 
togram) has systematically more numbers of galaxies in 
smaller b/a bins than the control galaxy sample (dotted 
histogram). For a better view of this over-density, we 
show the number ratios of the AGNs to the control galax- 
ies (Nagn/Nqal) in different b/a bins in Fig. [3] As one 
can see, the Nagn/Nqal decreases with the increasing 
of b/a systematically. The type 2 AGNs show both an 
excess of edge on objects and deficit of round objects. 
Since the galaxy sample have already been controlled to 
have the same physical properties and selection effects as 
the AGN hosts, the systemical bias in b/a distribution of 
the AGN hosts could only be stemmed from their pre- 
ferred (non-random) viewing angles. The bias towards 
the small axis ratios for the AGNs is basically consis- 
tent with the scenario that the nuclear obscurer is partly 
co-aligned with the stellar disk. 

3. THE SHAPE OF THE AGN HOSTS 

To quantify the preferred inclinations of the type 
2 AGN hosts, the distributions of their shape 
parameters (7,5) need to be known in advance. 

For the disk galaxies in general, it is found that 
their 7 distribution can be approximated by a nor- 
mal functio n whereas e follows a normal or log-normal 
distribution (jLambas et alJll992tlRvdeijl2004[ ). However, 
for our control sample of galaxies, the sample selection 
criteria is not well-defined and the goodness of fit is also 
bad when only single Gaussian distributions are assumed 
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Fig. 1. — The b/a, 7 and e distributions of the AGN hosts and 
control galaxy sample. Top panel: the b/a histograms of the AGN 
hosts(solid) and control galaxies(dotted). The dashed line shows 
the fitted histogram of the control galaxy sample from the non- 
negative least square regression. Bottom: The 7 and e distribu- 
tions reproduced from the non-negative linear regression, where 
the dotted lines show the contributions from different basis com- 
ponents (see Section 3). 

for 7 and e/lne. As an alternative, the control galaxy 
sample can be viewed as a combination of sub-samples 
of galaxies, where the 7 and e distributions both follow 
Gaussian distributions for each sub-sample. The 7 and 
e distributions of the control galaxy sample are then the 
sum of the Gaussian distributions of these sub-samples. 

We first assume 10 basis distributions for 7 and e each 
and then get 100 combinations of (7, e) distributions. 
Specifically, for ten 7 distributions, the dispersions are 
fixed to 0.02, while the mean values are assigned from 
0.1 to 0.5 with a step of 0.04. For e, the peaks of the 
Gaussian distribution are fixed to and the dispersions 
are taken values from 0.02 to 0.2 with a step of 0.02. Each 
of the (7, e) combination is then viewed as a sub-sample. 
Using Monte- Carlo simulations of random viewing an- 
gles, we obtain the b/a distributions of each sub-sample. 
Finally, we use the non-negative least square linear re- 
gression technic to fit the observed b/a distribution of the 
control galaxies to get the fraction of each sub-sample. 
The final fitted b/a is shown as the dashed line in Fig. m 
As we can see, the fitting is very well, which is consistent 
with the observed value inside Possion error in most of 
the b/a bins. The resulted (7, e) distributions of the con- 
trol galaxy sample are shown in the bottom two panels 
of Fig. [TJ where the dotted lines in each panel show the 
contributions from different sub-samples. 

With (7, e) distributions known for the control galax- 
ies, so that also for the AGN hosts, we move to model 
the viewing angles of AGN hosts. 

4. THE ORIENTATION OF THE NUCLEAR 
OBSCURER 

The biased axis ratio distribution (towards small val- 
ues) of the AGN hosts implies that either the orienta- 



orientation of nuclear obscurer 
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Fig. 2. — The cartoon of the tilt angle (A<rg) between the torus 
plane (inclination: iy) and stellar disk (inclination: ig) of a type 
2 AGN (90° -% T < e T ). 

tion of the molecular torus is correlated with the stellar 
disk or there is an extra dust layer aligned with stellar 
disk. We denote these two different obscuring scenar- 
ios as the 'aligned torus model' and 'stellar dust model' 
respectively below. 

We parameterize the inclination of the dust torus as it, 
where «t (0° < %t < 90°) is the angle between the normal 
of the torus plane and the line-of-sight. Correspondingly, 
the inclination of the stellar disk is denoted as is. The 
tilt angle between the torus plane and stellar disk is then 
denoted as Ats- 

4.1. the obscurer: aligned torus model 

In the aligned torus model, we assume that the obscur- 
ing torus has an opening angle 2©t- Obviously, a type 
2 AGN satisfies the criteria 90° — %t < ©t- A cartoon of 
such a torus model is shown in Fig. [2] The ©t and the 
tilt angle Ats are then the two free parameters in this 
model. 

©t could be constrained from the type 2 AGN fraction 
/2 since = sm(6r). The studies on the local Seyfert 
galaxies and the statistics of active galaxies in SDSS have 
shown that /2 is in th e range o f 70-9 percent for the 
low luminosity AGNs ()Ho et all 119971 : iSimpsonI 120051 ) , 
corresponding to 45° < ©t < 65°. Thus, we take 45° 
and 65° as the low and up limits of Ot- Inside this 
range, Ot is set as a free parameter. 

For tilt angle Ats, it is a reasonable assumption that 
the most probable orientation of the torus plane is to be 
the same as the stellar disk, i.e. P m ax(is) = P(Ats = 0) 
, but with significant scatters on Ats- On the other 
hand, cos(Ats) would follow a uniform distribution be- 
tween and 1 if the orientations of the torus plane and 
the stellar disk were independent. Thus, we further 
assume that cos(Ats) follows a Gaussian distribution 
which peaks at Ats = and has a scatter cos(ATs,m)- 



Under this assumption, Ats,™ could be viewed as an ef- 
fective tilt angle, within which 68 percent of the Ats are 
smaller than. 

With ©t and Ats,™ assigned, we make Monte-Carlo 
simulations to predict the number ratios Nqal/Ngal m 
different b/a bins. In specific, we first generate a large 
Monte-Carlo sample of galaxies with random viewing an- 
gles and (7, e) distributions as shown in Fig. [TJ For each 
galaxy with viewing angle (is,it4>s,i)i we get the incli- 
nation of the torus plane iT,i from the tilt angle Ats,i 
by assuming a random 4>T,i, where <px,i and <ps,i are the 
position angles of the torus plane and stellar disk respec- 
tively. After that, we get a modelling type 2 AGN sample 
by applying the criteria 90° — its < ©t for all Monte- 
Carlo galaxies. The other galaxies with 90° — %t.i > ©t 
are classified as type 1 AGNs. 

In the observed type 2 AGN sample, there are con- 
taminations of intermediate type AGNs, which show 
weak broad wings on the Ha lines but are classi- 
fied as 'galaxies' by the SDSS spectroscopic pipeline. 
The fraction of these intermediate type contamina- 
tions is estim ated to be about 8 per c ent from eyeba ll 
classifications(Kauffmann et al. 2003|; ISimpsonI [2005). 
To mimic this intermediate type AGN contaminations 
in our model as well, we also put some type 1 AGNs into 
our type 2 AGN sample in random as contaminations, 
up to the maximum of 8 percent of the total sample. 

For each Monte-Carlo simulation, we randomly se- 
lect 50,000 objects from the resulted library of type 2 
AGNs(including type 1 contamination). For each param- 
eter set of ©t and A T ^ m , we run the simulation 30 times 
to account for the Possion fluctuation. We normalize the 
b/a histograms of the model galaxies to be the same as 
the model AGNs and calculate Nagn/Ngal m b/a bins 
as that done for observations. The mean values of the 30 
simulations are taken as the model output for each pa- 
rameter set. Finally, we use the least square routine to 
search the best model parameters of ©t and Ats,™ by 
fitting the observed Nagn /Nqal in different b/a bins. 
The best model fitting is shown as the dashed line in 
Fig. [3j which has the minimum Xmin ~ 9, comparable 
to the number of data points(9 b/a bins), indicating an 
excellent goodness of fit. 

The best model parameters are ©t = 58° and 
Ars.m = 37° respectively, ©t = 58° corresponds to a 
type 2 AGN fraction of [2 = 0.85, in excellent agreement 
with the res ults quoted in literature for t he low lumi- 
nosity AGNs l|SimpsoiJl200HlLu et al.l[2f)Tfl . For the tilt 
angle between stellar disk and torus plane, the effective 
tilt angle Ats,™ = 37° corresponds to a mean tilt angle 
Ats ~ 30°. We emphasis here that our assumption of 
the Gaussian distribution of cos Ats is in a arbitrary way, 
however, the mean tilt angle (Ats ~ 30°) is quite robust 
and independent of the shape of the Ats distribution. 

We show the 68 (A x 2 = 2.30) and 90 (A x 2 = 4.61) 
percent of the confidence levels of our model parameters 
in the sub-panel of Fig. [3] As we can see, the parameter 
©t is degenerated with Ats,™ to a certain extent. The 
effect of increasing the torus opening angle ©t is similar 
to the increasing of the effective tilt angle Ats,™, which 
both result in a more randomized viewing angles of the 
disks of AGN hosts. Therefore, an independent measure- 
ment of ©t from the type 2 AGN fraction is crucial to 
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Fig. 3. — The ratio of the numbers of the AGNs to the control 
galaxies in different b/a bins. The solid line with error-bars shows 
the observational results from the SDSS. The dashed curve shows 
the best model fitting from the aligned torus model(see Section 
14. Il l while the dotted curve show the model prediction from the 
stellar dust model(see Section |4,2|I , The dot-dashed horizontal line 
is constant 1. The sub-panel shows the confidence level of the two 
model parameters, 9y and A^g m for the aligned torus model. 



constrain the tilt angle Ats in our model. 

4.2. the obscurer: stellar dust model 

In this section, we test the scenario that the over- 
density of the edge-on galaxies of type 2 AGNs is caused 
by an extra dust layer aligned with the stellar disk. 
Therefore, we no more assume that the nuclear obscurer 
is aligned with the stellar disk as that in the aligned torus 
modelfSection ETTT) . but consider it being randomly ori- 
ented. Besides the random orientated torus, the dust 
aligned with the stellar disk will obscure the nucleus 
when it is viewed as edge-on. 

Using the denotations as in aligned torus model, a ran- 
domly oriented nuclear torus means that the orientations 
of the stellar disk and torus plane are independent of 
each other, i.e. cos(Ats) follows a random distribution 
between and 1. For its opening angle, we set Ot = 45°, 
the lower limit of the aligned torus model, since we have 
an extra component of type 2 AGNs obscured by the 
host-aligned dust layer. For host-aligned obscurer, we 
assume its height is the same as the stellar disk. Thus, 
for a stellar disk with height 7, its nucleus will be ob- 
scured when its viewing angles satisfies tan(is)< 7. 

We then combine these two types of obscurers and get 
the final type 2 AGN sample. The other model set- 
tings and procedures are the same as in aligned tours 



model. We show the model prediction as the dashed 
line in Fig. [3l As we can see, this extra stellar ob- 
scurer model indeed predicts over-density of the low axis 
ratio galaxies for the AGN sample. However, this over- 
density mostly takes place at b/a < 0.4. For the others 
with b/a > 0.4, the predicted Nagn/Nqal is almost a 
constant, whereas the observations show continuous de- 
creasing of type 2 AGNs in more round galaxies. This is 
because that most of the stellar disks having 7 < 0.4(see 
Fig. [T|). Therefore, this model prediction shows signif- 
icantly poorer fitting(x 2 w 36, outside the 99.9% con- 
fidence interval of the x 2 distribution with 8 freedoms) 
than the aligned torus model. 

However, we would mention that there is no tune- 
able parameter in this stellar dust model, whereas the 
aligned torus model has two free parameters. Setting 
8t to be a free parameter in the stellar dust model does 
not help, because the torus was assumed to be orien- 
tated randomly. However, assuming that the stellar dust 
layer following the same shape of the stars might be too 
simplistic. A more sophisticated treatment of the dis- 
tribution of the s t ellar dust might improve the fit (e.g. 
IMaiolino fc Riekel (|1995l )1. but is out of the range of this 
study. 

5. CONCLUSION AND DISCUSSION 

In this study, we select a sample of type 2 AGNs hosted 
by spiral galaxies from the DR7 of the SDSS. We build a 
control sample of non-active galaxies to the AGN sample 
by matching their observational and physical properties 
in detail. By comparing with the control sample, we find 
that the AGN hosts are systematically biased to edge- 
on views. By modelling the height and ellipticity of the 
disk galaxies with a non-parametric technic, we further 
find that this systematical bias can not be quantified by 
an extra dust layer aligned with stellar disk but can be 
nicely reproduced by assuming an average tilt angle of 
~ 30 degree between the torus plane and stellar disk. 

Some of the assumptions in our model are simplified, 
e.g. the constant type 2 AGN fraction. There are stud- 
ies suggesting that the opening angle of the torus is 
correlated with the central AGN luminosity. For high 
luminosity AGNs, the distance of the torus from the 
nucleus is further ('receding torus model') so that the 
torus o pening angel and the type 2 AGN fraction are 
smaller (|Lawrenc e f[l99lTl . Most of our AGNs are low lu- 
minosity objects(L[oni] < 10 8 L©), and the type 2 AGN 
fraction has been shown to be roughly a constan t in this 
luminosity region (|Simpsodl2005t ILu et al.|[2010h . A sep- 
aration of our sample into two sub-samples according to 
-^[oiii] also would not show any significant differences but 
reduce the goodness of the statistics. 

In the unified AGN diagram, the torus lies on the same 
plane of the accretion disk. In this case, the mean tilt 
angle between the torus plane and the stellar disk con- 
strained from our model could also be understood as the 
mean tilt between the accretion disk and stellar disk. 
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